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ABSTRACT 

We present the time variability properties of a sample of six blazars, AO 0235+164, 3C 273, 3C 
279, PKS 1510-089, PKS 2155-304, and 3C 454.3, at optical-IR frequencies as well as 7-ray energies. 
These observations were carried out as a part of the Yale/SMARTS program during 2008-2010 that has 
followed the variations in emission of the bright Fermi-LAT- monitored blazars in the southern sky with 
closely-spaced observations at BVRJK bands. We find that the optical-near IR variability properties 
are remarkably similar to those at the 7-ray energies. The discrete auto-correlation functions of 
the variability of these six blazars at optical-IR and 7-ray energies do not show any periodicity or 
characteristic timescale. The power spectral density (PSD) functions of the R-band variability of all 
six blazars are fit well by simple power-law functions with negative slope such that there is higher 
amplitude variability on longer timescales. No clear break is identified in the PSD of any of the 
sources. The average slope of the PSD of R-band variability of these blazars is similar to what was 
found by the Fermi team for the 7-ray variability of a larger sample of bright blazars. This is consistent 
with leptonic models where the optical-IR and 7-ray emission is generated by the same population 
of electrons through synchrotron and inverse-Compton processes, respectively. The prominent flares 
present in the optical-IR as well as the 7-ray light curves of these blazars are predominantly symmetric, 
i.e., have similar rise and decay timescales, indicating that the long-term variability is dominated by the 
crossing time of radiation or a disturbance through the emission region rather than by the acceleration 
or energy-loss timescales of the radiating electrons. For the blazar 3C 454.3, which has the highest- 
quality light curves, the total energy output, the ratio of 7-ray to optical energy output, and the 
7-ray vs. optical flux relation differ in the six individual flares observed between 2009 August and 
December. The results are consistent with the location of a large 7-ray outburst in 3C 454.3 during 
2009 December being in the jet at ~18 pc from the central engine. This poses strong constraints on 
the models of high energy emission in the jets of blazars. 

Subject headings: black hole physics — radiation mechanisms: non-thermal — galaxies: active - 
galaxies: jets — (galaxies:) quasars: general — (galaxies:) quasars: individual: 3C 
454.3 



1. INTRODUCTION 

Blazars are a sub-class of active galactic nuclei (AGN) 
characterized by a promine nt jet pointing within a few 
degrees of our line of sight (|Urrv fc Padovanilll995D . In 
many cases, the jets are luminous over a wide range of 
wavelengths from radio to 7-rays. Due to the proximity 
of the jet axis to the line of sight, emission from the jet is 
relativistically beamed and hence amplified by an order 
of magnitude or more for many blazars. The observed 
timescales are also shorter than that in the rest frame 
of the jet plasma. Most of the observed radio to optical 
(and in some cases X-ray) emission fro m the blazars is 
due to synchrotron radiat i on in the jet (|Bregman et al 



1981 1 : lUrrv k M ushotzkvl 119821: llmpev fc Neugenbauer 
T988tlMarscherlll998D . X -rays and 7-rays may be due to 



inverse Compton scattering by the same energetic elec- 
trons radiating synchrotron emission (the so called lep- 
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tonic models; e.g., [Bottchcr 2007J) or due to synchrotron 
radiation by protons co-accelerated with the electrons in 
the jet, interactions of these highly relativistic protons 
with external radiation fields, or p roton-induced parti- 
cle cascades (hadronic m odels; e.g. IMucke fc Protherog 
l2lMlMuck"e et al.l[200l . 

In the leptonic model, the source of the seed photons 
that are being scattered may be the synchrotron photons 
generated within the jet, in which case i t is termed syn- 
chrotron self-Compton (SSC) process (IMaraschi et al.l 
Il99l [Chiang fc Bottcherl [20021: lArbeiter et all 120051) . or 
from outside the jet (radiation from broad emission 
line region or BLR, accretion disk, or dusty torus), 
term e d external Compton or EC process (Sikor aTet all 
199l[Coppi fc Aharonianlll999t iBlazeiowski et al.ll2000t 



Dermer et alll2009D . 



Time variability across multiple wavebands is a defin- 
ing property of blazars and has been used to probe 
the location and physical processes rel ated to the emis- 
sion at very fine resoluti ons (e. g.. iChatteriee et all 

20081 iMarscher et al.l [200l [20Toi: Uorstad et all 120101: 



Agudo et al J [2Q11J). Most of the observed radiant en- 



ergy produced in many blazars peaks in the 7-ray part 
of the spectrum. But until recently, long-term and well 
sampled 7-ray light curves of blazars were not avail- 
able. Expansion of such monitoring to a wide range of 
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7-ray energies (20 MeV to 300 GeV) is now occurring 
through the Fermi Ga mma-Ray Space Te lescope which 
was launched in 2008 (|Abdo et alj|2010b[ ). The Fermi 
bright so urce list from the fir st 11 months contains more 
than 600 (lAbdo et al.ll2010al ) and the Fermi 2-yr catalog 
(jAbdo et al.H2011| ) contains more than 1000 AGNs most 
of which are blazars. 

Investigating the time variability of the data from 
Fermi and supporting multi-frequency monitoring with 
statistically robust techniques is one of the most promis- 
ing ways to understanding how and where the multi- 
wavelength emission is produced in the jets of blazars. 
For example, in the leptonic models, the low and high- 
energy variability should be correlated while in some ver- 
sions of the hadronic models, this correlation may be 
weak or absent. For this reason, multi- frequency moni- 
toring programs are crucially important. 

Here we report on the Yale/SMARTS optical-near IR 
monitoring program^ during 2008-2010 which has fol- 
lowed the variations in emission of the Fermi-LAT mon- 
itored blazars in the southern sky w ith closely-spaced 
observations (see lBonning et al~ll2012L for details of data 
acquisition and calibration). Six blazars were detected 
for a sufficiently long interval by the Fermi-LAT during 
2008-10: AO 0235+164, 3C 273, 3C 279, PKS 1510-089, 
PKS 2155-304, and 3C 454.3. In order to quantify their 
variability properties, we calculate the power spectral 
density and the auto-correlation function, and carry out 
flare decomposition analysis of the SMARTS light curves 
of these six sources. We then compare the results with 
similar properties of the 7-ray variability as discussed by 
lAbdo et all (|2010bD who analyzed 11-month-long GeV 
light curves of a sample of bright blazars from the Fermi 
3-month source catalog. Possible models of emission may 
be distinguished by comparing the characteristics of the 
flux variations at lower (e.g. optical-IR) and higher (e.g. 
7-ray) energies. 

The auto-correlation function (A CF) can provide in - 
sights into the nature of variability. lAbdo et al.l ()2010bl ) 
have shown that the ACFs of the 7-ray variability of a 
sample of bright blazars do not show any hint of quasi- 
periodicity or characteristic timescale. Here we calcu- 
late the ACFs of those six blazars to search for such a 
timescale in optical or near IR bands. Comparison of the 
ACF of the optical-IR variability with that of the 7-ray 
variability can be used to understand the relation be- 
tween the mechanism and location of emission in these 
wave bands. For example, in the leptonic models, the 
optical-IR and 7-ray emission is generated by the same 
population of electrons. Hence, it is expected that the 
width and shape of the ACF will be very similar at both 
energies. However, some hadronic models are also able 
to reproduce this behavior and hence we cannot exclude 
those models on this basis. In this paper, we calculate 
the ACFs of the 7-ray variability of those blazars as well 
to facilitate the comparison of the 7-ray ACFs to those 
at optical-IR wave bands. 

Power spectral density (PSD) analysis is a common 
technique to characterize time variability. The PSD 
corresponds to the power in the variabil it y of em ission 
as a function of timescale. lAbdo et al.l (|2010bD show 
that the PSDs of the 7-ray light curves of a sample of 

5 http:/ /www. astro. yale.edu/smarts/fcrmi/ 



bright blazars are well-described by simple power laws, 
corresponding to "red noise" . Red noise is defined as 
uncorrelated fluctuations where power density decreases 
with increasing frequency. In the case of an astronomical 
time series this translates to having larger amplitude 
variations at longer timescales. X-ray PSDs of Seyfert 
galaxies and black hole X-ray binaries (BHXRBs) can be 
fit by a piece-wise power law, with a slope that steepens 
above a "break-frequency" . The break-frequency is 
related to a characteristic timescale (the origin of 
which is not well-understood) that is proportional to 
the mass of the central black hole (No wak et al.l Il99£ 



Uttlev. McHardv. fc Papadakisl l200a 



2004t iMarkowitz et al. l 120031: 



.Pounds et alj 

Edelson fe Nandral 119991 : iChatteriee et all I200S 



We calculate the PSD of the optical variation in these 
blazars in order to search for the existence of such a 
characteristic timescale. Similar to the ACFs, if the 
optical and 7-ray emission is radiated by the same 
electrons, the shape and slope of the PSD functions at 
these two energy-bands are expected to be similar. 

Where well-sampled simultaneous light curves exist, 
we can compare the properties of the individual 7-ray 
and optical-IR outbursts to investigate the nature of the 
emission mechanisms. The total energy output of an 
optical-IR flare can also be compared to that of the cor- 
responding 7-ray flare. To achieve this, we decompose 
the light curves i nto individual ( somet imes overlapping) 
flares, similar to lValtaoia et al.1 (119991), each with expo- 
nential rise and decay. lAbdo et al.1 (|2010b[ ) have also 
analyzed segments of Fermi light curves of 10 sources 
to investigate the properties of individual flares. They 
found that most of the outbursts are symmetric in na- 
ture. In this paper, we use a slightly different algorithm 
to perform a similar analysis of the SMARTS R-band 
and Fermi light curves of the six blazars in our sample. 

There are two main classes of blazars, namely, flat 
spectrum radio quasars (FSRQs) and BL Lacertae type 
objects (BL Lacs). One of the main differences in the 
observable properties of these two classes is that the BL 
Lacs have much weaker or no observable emission lines. 
This is probably due to different parent populations of 
this two classes, i.e., FSRQs are linke d with FR II and 
BL L acs are with FR I radio galaxies (jUrrv fe PadovanH 
119951 ). Although the bright jet continuum emission com- 
plicates our view of the unbeamed thermal components 
(jGeorganopoulos fe Marschen |1998|) . and it has caused 
some deba te about the meaning of the BL Lac-FSRQ 
difference (IVermeulen et al.l 119951: ICorbett et all 12000 : 
IRaiteri et al.l I2007t) . it is now apparent that the ma- 
jority of BL Lac objects have intrinsically less luminous 
thermal emission, pointing at some fundamental dif- 
ference (lUrrv fe PadovanH Il995t iFossati fe Meveil 120101: 
IGhisellini et al.ll201lD . A change in the properties of their 
central regions associated with a critical value of the mass 
accretion rate (i n Eddington units) has been proposed 
(|Ghisellini et al.ll2009D . and there are indications that it 
might hold for radio-g alaxies, too (IGhisellini et al.ll2010l : 
iChiaberge et al.ll2002D . 

In addition to the FSRQ/BL Lac divide, the fre- 
quency of the synchrotron peak (in vF u vs. v diagrams) 
has emerged as one of the most important observa- 
tional distinctions, leading to the classification of blazars 
based on its value, namely, Low/Intermediate/High Syn- 
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Table 1 

The half width at half maximum (HWHM) 
of the autocorrelation functions of the 
blazars at three wave bands. 



HWHM (Days) 



Object 


B-Band 


J-Band 


7-ray 


AO 0235+164 


24 


24 


26 


3C 273 


14 


13 


8 


3C 279 


41 


28 


8 


PKS 1510-089 


8 


12 


5 


PKS 2155-304 


27 


31 


11 


3C 454.3 


9 


9 


9 



chrotr on Peak objects (LSP/ISP/HSP). iFossati et al.l 
(1998) showed that blazar SEDs seem to change sys- 
tematically with luminosity: the most powerful ob- 
jects are LSP, while HSP SEDs are associated with 
relatively weak sources. This result has been sup- 
ported by studies of high-z blazars, low power objects 
and r ec ently of 7-ray selected samples dFabian et al.l 
2001albt IC ostaman te et aLl|2001i iGhisellini fc Tavecchiol 
2009HSambruna et al.ll2010h . but more recent work sug- 
gests that the hypothesized relationship between SED 
peak and lu minosity may be m ore complex than a sim- 
ple sequence (|Mever et al.ll20lilh 

Among the six blazars in our sample, PKS 2155-304 
is a high-synchrotron-peaked BL Lac type object and 
the rest are FSRQs as well as low-synchrotron-peaked 
objects. Here we test if any significant difference exists 
between the multi-wave band variability properties of the 
FSRQs in our sample and that of PKS 2155-304. 

In §2 we calculate the auto-correlation functions of the 
optical, infrared, and 7-ray light curves of our sample 
of six blazars while in §3 we present the power spectral 
analysis, its results and implications. In §4 we determine 
the shortest variability timescales present in our data. 
We model the flares in optical and 7-ray light curves 
and discuss the properties of the flares in §5 while in 
§6 we compare the time variability properties of FSRQs 
and BL Lacs. In §7 we decompose a segment each of 
contemporaneous Fermi and SMARTS light curves of the 
blazar 3C 454.3 to investigate the nature of the outbursts 
and the corresponding emission processes. §8 presents 
the summary and conclusions. The 7-ray light curves 
analyzed in this work are at 0.1-300 GeV energy band 
and are taken from the Fermi LAT table of monitored 
sources provided by the Fermi team. 

2. AUTO-CORRELATION FUNCTION 



The discrete cross-correlation function (DCCF) is a 
method to calculate t he cross-correlation of unevenly 
sampled discrete data (jEdelson &; Kr olik 1988]). We use 
this to calculate the discrete auto-correlation function 
(DACF) of the variability of the blazars in our sample at 
7-ray energies, as well as B and J-band frequencies. The 
width of the ACF may be related to a characteristic size- 
scale of the corresponding emission region, while equally 
spaced and repeated peaks or drops in the function shape 
can point out characteristic timescales and provide hints 
of possible quasi-periodicity. 

The optical-IR light curves obtained by SMARTS have 
seasonal gaps (~3-4 months) due to the source going be- 



hind the Sun. These large gaps distort the ACF. To avoid 
this, for each of the six blazars, we select the longest seg- 
ment without such large gaps and use them to calculate 
the ACFs. The length of such segments were 200—250 
days for all six blazars in our sample. This enables us 
to avoid distortion due to large data-gaps while keep- 
ing the resultant ACF a significant representative of the 
nature of the entire light curve. The average sampling 
rate of the SMARTS light curves was one data point 
every 2-3 days. The light curves from Fermi do not con- 
tain the Sun-gaps but the blazar may not be consistently 
detected during the selected segments of the optical-IR 
light curves. Therefore, while making the selection of the 
optical-IR segments, we chose those intervals when the 
blazar was detected by Fermi at least once in 2-3 days 
on average so that the 7-ray sampling is similar or bet- 
ter than that in the optical-IR bands. We chose the B 
and the J-band light curves for comparison with the 7- 
rays since those are the bands with the widest separation 
among our data while the K-band data have a sampling 
frequency less than once in 2-3 days in some cases. 

The ACF of a light curve depends on the temporal sam- 
pling as well as the intrinsic variability of the source. The 
sampling frequency in the J-band is better than that in 
the B-band. Hence, we resampled the 7-ray and J-band 
light curves with that in the B-band. To carry out the 
resampling, we remove data points from the 7-ray and 
J-band light curve segments for which there are no cor- 
responding B-band data point within ±0.5 days. This 
was done to ensure that the comparison of the ACFs is 
physically meaningful. In 3 out of 6 blazars in our sam- 
ple, namely AO 0235+164, PKS 1510-089, and 3C 454.3, 
the B-band as well as the resampled 7-ray and J-band 
light curves had a sampling of at least once in 3 days. 
The ACFs of those 3 sources at all three wave bands are 
shown in Figure [TJ For the ACF calculation, we use a 
bin-width of 5 days. We list the half-width at half maxi- 
mum (HWHM) of the central peaks of the ACFs in Table 
[TJ The ACFs in all bands are very similar in shape for 
AO 0235+164 and 3C 454.3. In PKS 1510-089, the 7- 
ray ACF is slightly narrower in shape but the HWHM 
is very similar to the that in B and J-bands. Since the 
temporal sampling in all 3 bands are the same, this is a 
comparison of the intrinsic variability properties of these 
three blazars. This shows that the characteristic size- 
scales of the corresponding emission regions in all three 
bands are similar which, in turn, may imply that the 7- 
ray and optical-IR emission in these blazars is generated 
by electrons of similar energies. Even initially homoge- 
neous electron acceleration region will develop energy- 
dependent stratification due to the energy dependence 
of radiative loss timescales, i.e., lower energy electrons 
naturally diffuse into a larger volume. Similarly, a higher 
energy population of electrons will be in a smaller volume 
causing the ACF width and characteristic timescales of 
variability to be smaller. 

Figure [U presents the ACFs of the other 3 blazars, 
namely, 3C 273, 3C 279 and PKS 2155-304, for which 
resampled data were not used due to poor sampling. In 
3C 273, the ACFs at different bands are similar while in 
3C 279, the 7-ray ACF is significantly narrower. This 
may be an effect of the temporal sampling and can not 
be confirmed with the present data. The peak of the 7- 
ray ACF of the blazar 2155-304 is much smaller than 
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Figure 1. Discrete auto-correlation functions of the 7-ray (blue dotted line), B-band (red solid line and filled circles with error bars), and 
J-band (green dashed line) light curves of AO 0235+164, PKS 1510-089, and 3C 454.3. Uncertainties shown by the red data points are 
characteristic of other bands as well. In these three sources, the temporal sampling in all 3 bands is the same, hence this is a comparison 
of the intrinsic variability properties of these blazars. 
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Figure 2. Discrete auto-correlation functions of the 7-ray (blue dotted line), S-band (red solid line and filled circles with error bars), 
and J-band (green dashed line) light curves of 3C 273, 3C 279, and PKS 2155-304. Because of the scarcity of data, sampling differed at 
optical/IR and 7-ray wavelengths. Uncertainties shown by the red data points are characteristic of other bands as well. 



the rest. This is because, the 7-ray variation of this 
source was dominated by short-term flares uncorrelated 
to each other. This behavior is also reflected in the power 
spectrum of its variation as discussed in §3. We do not 
find any significant peak in any of the ACFs other than 
the central peak. That implies the absence of significant 
quasi-periodic behavior in the variability of these blazars 
in 7-ray and optical-IR frequencies. 

3. POWER SPECTRAL ANALYSIS 

We use a varia nt of the Power Spectrum Response 
method fPSRESP: IUttlev. McHardv. fc Papadakisl2002t ) 
to determine the intrinsic PSD of the optical light 
curves. Our realizati on of PSRESP is described in 
iChatteriee et a l. (2008) in details. Here we briefly sum- 
marize the method. 

We start with an underlying model for the PSD func- 
tion, such as a simple power-law. Then we simulate N 
light curves starting fro m this underlying model using 
the technique described in lTimmer fc Koenigl ()1995l ). We 
use N = 100. We resample these simulated light curves 
with the sampling window of the observed light curve. 
We also add other noise to the simulated light curves 
which affect the real observations due to its finite length 
("red noise leak") and discontinuous sampling ("alias- 
ing"). Then we calculate the PSD of the real light curve 
(denoted by PSD Q b s ) and that of each of the simulated 
light curves (denoted by PSD s j mi ;, i=l, N). We calcu- 



late the approximate value of the power due to Poisson 
noise using the observational uncertainties and add that 
power to PSD S i m ,j. We finally compare PSD Q b s with the 
mean of the PSD s j mi ; (denoted by PSD S i m ) weighted by 
the standard deviation in the distribution PSD S i mj ; to 
determine the goodness of fit of the underlying model 
that we assumed at the beginning. The goodness of fit is 
quantified by the "success fraction" F succ defined as the 
fraction of simulated light curves that successfully rep- 
resent the observed light curve. We do this for a range 
of underlying models covering a large parameter space 
and determine the underlying model that provides the 
best-fit, i.e., the highest F succ . The observed PSD suf- 
fers from the distorting effects of the finite length ("red 
noise leak") and discontinuous sampling ("aliasing") of 
the light curves as well as power generated by the tem- 
poral sampling. These are accounted for by the method 
that we use. 

Similar to the ACF calculation, we select the longest 
segments of SMARTS R-band light curves of the six 
blazars without seasonal gaps (~3-4 months). This is 
done in order to minimize the distortion of the power 
spectrum and to cover the largest possible range of 
timescales. The length of such segments were 200—250 
days for all six blazars in our sample and the average 
sampling rate was one data point every 2-3 days. We 
chose the R-band for this calculation because among the 
optical bands it has the best sampling for all objects. 
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Figure 3. Power spectral density of the i?-band light curves of the six blazars in our sample (calculated using the PSRESP method). 
The PSD of the observed data is given by the solid jagged line, while the underlying power-law model is given by the dashed straight line. 
Filled circles with error bars correspond to the mean value of the PSD simulated from the underlying power-law model. The error bars are 
one standard deviation of the distribution of simulated PSDs. The slopes of the best-fit models are listed in Table [2] All PSDs are well- fit 
by simple power-law models with no break. 
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Figure 4. The data points show the average power spectral den- 
sity of the six blazars in our sample. The error bars are one stan- 
dard deviation about the mean. The average slope —1.6 ± 0.2 is 
denoted by the dotted line. 

Figure [3] presents the PSDs of the six blazars and 
the corresponding power-law fits. The solid jagged line 
shown in Figure [3] is PSD Q b s and the filled circles with 
error bars are the mean (PSD s j m ) and standard devia- 
tion of the distribution PSD s j m .i. The resampling and 
addition of the noise described above add extra power 
to PSD s i m causing the filled circles with error bars to be 
slightly higher than the underlying model (dashed line) 
at certain frequencies. 



The optical R-band PSD of all blazars in our sample 
show red noise behavior, i.e., there is higher amplitude 
variability on longer than on shorter timescales. We list 
the best-fit slope and success fraction of the PSD of all 
blazars in Table The rejection confidence, equal to 
one minus the success fraction, is much less than 0.9 in 
all cases. This implies that a simple power-law model 
provides an acceptable fit to all the PSDs. We also fit 
a broken power-law model to all the PSDs, setting the 
low- frequency slope at —1.0 and allowing the break fre- 
quency and the slope above the break over a wide range 
of parameters (10~ 7 to 10~ 6 Hz and —1.0 to —2.5, re- 
spectively) while calculating the success fractions. This 
gives lower success fractions than the simple power-law 
model across the entire parameter space. Hence, there 
is no significant break detected in the PSD of any of the 
sources. We determine the optical PSD up to the high- 
est variational frequency that can be achieved with the 
existing data. A better constraint on the existence of 
a break in the optical PSD might be achieved with a 
broader range of sampled frequencies. 

Among the blazars in our sample, the i?-band vari- 
ability amplitude of 3C 273 was much smaller than the 
others during the time interval considered here. This is 
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Table 2 

Slope and Success Fraction of the Best-Fit Power Spectral 
Models. 



Object 


Slope 


succ 


AO 0235+164 


-1 6+ U ' 3 


0.74 


3C 273 


o o+0.2 
z - J -0.5 


0.63 


3C 279 


9 O+0.5 
z,o -0.2 


0.99 


PKS 1510-089 


-0 6 +0 ' 5 

U.D_q 2 


0.87 


PKS 2155-304 


9 9 + 0.2 
z ' z -0.4 


0.93 


3C 454.3 


1 B+0.3 


0.81 



evident in the PSD which shows that the power of vari- 
ability in 3C 273 is about two orders of magnitude less 
than the average of the other blazars at all timescales. 
The PSD of the blazar PKS 1510-089 has a significantly 
flatter slope than the other blazars in our sample indicat- 
ing larger amplitude of short-timescale variations. This 
trend is evident in the light curve of this blazar which is 
dominated by very large-amplitude flares of width a few 
days. 

The average power spectral density of the six blazars 
in our sample is shown in Figure 01 The data points 
and the error bars are the mean and standard devia- 
tion of the power of all six blazars in each logarithmic 
bin. The best fit slope of the average PSD is — 1.6 ± 0.3, 
where the uncertainty i s one s tandard deviation about 
the mean. lAbdo et all (|2010bD found that the average 
slope of the best-fit power-law for a sample of 9 bright 
FSRQs is -1.4 ± 0.1 and that of 6 bright BL Lacs is 
— 1.7 ±0.3, similar to the slope we found. This is consis- 
tent with the model where the emission at 7-ray energies 
and optical frequencies are both generated by electrons 
through synchrotron and inverse- Compton processes, re- 
spectively, so that the nature of variability at these two 
bands is similar. The uncertainty in the individual as 
well as the average PSD will decrease when more data 
are included and may reveal the existence of small dif- 
ference in the amplitude of power between the 7-ray and 
optical-IR variability at a given frequency. For example, 
the optical PSD of the BL Lac object PKS 2155-304 is 
quite steep (slope ~ — 2.2) while the 7-ra y PSD is much 
flatter (see Fig. 17 of lAbdo et aHl2010b| ). This implies 
that the short-term (few to 10 days) variability in this 
blazar is more prominent at 7-ray energies than at opti- 
cal wave bands. Possible explanation for this is discussed 
in §6. 

4. SHORTEST TIMESCALES OF VARIABILITY 

We calculate the doubling timescales (Td ou b) for each 
object in all bands. Tdoub is defined as the minimum 
time interval over which the flux increases by a factor of 
2. We scan the light curves of all objects in our sample in 
each band to find time intervals over which the flux has 
doubled and then select the minimum of those time in- 
tervals as the doubling timescale for the respective band. 
Similarly, we also determine the minimum timescale over 
which the flux has decreased by a factor of 2 and denote 
it by Thaif- To avoid spurious results, we select Td ou b 
and Thaif only in those cases where increase or decrease 
of flux at that timescale has happened 3 or more times 
during the period we scanned and the general trend is 
present at all bands. 



Table [3] lists Td ou b and 7h a if for all objects in our sam- 
ple in all available bands. It can be seen from the table 
that the timescales over which the flux increases or de- 
creases by a factor of 2 is similar in a given band for an 
object. This broad similarity implies that the flares are 
symmetric and hence crossing time dominated. However, 
the value of Thaif is larger than Tdoub in 3C 454.3. This 
blazar was not observed by SMARTS during its decay 
from the large outburst in 2009 December due to Sun- 
gap. The faster decay and hence a smaller value of Thaif 
could be detected during that time. 

No object shows a significant difference in those 
timescales among the optical-near IR bands. However, 
in 3C 454.3, the doubling timescale is smaller at longer 
wavelength optical-IR bands. This may be due to the 
disk contribution in the bluer bands (e.g.. iBonnoli et all 
120111) which dilutes the jet variability and hence increases 
the doubling timescale. This is supported by observation 
of 3C 454.3 at B, V, R, and J bands during 2009 Decem- 
ber as shown in Figure [5] The B-band flux changes by a 
factor of only ~1.5 over 2-4 days while the corresponding 
change in J-band is by a factor of ~6. 

Tdoub and Thaif of all 6 objects in the 7-ray band is ~1.0 
day. These are all upper limits since the timescales are 
determined from light curves provided by the Fermi team 
which are binned at 1-day intervals. Those timescales 
for AO 0235+164 and PKS 1510-089 are ~2 days in 
all optical-IR bands, consistent with those in the 7-ray 
band. These are also upper limits since we do not have 
data points between those 2 days to check if the actual 
timescales are significantly less than that. In the case of 
3C 273, the flux value does not change by a factor of 2 
over any timescale during the two yr of monitoring. In 
3C 279, those timescales at all optical-IR bands is ~4-5 
days. This is not an upper limit since we do have data 
points during those 4 days and the flux did not change 
by factor of 2 in less than 4 days. The slower variability 
indicated by the longer doubling timescales in optical-IR 
bands than that in the 7-rays is consistent with the nar- 
rower auto-correlation function of 3C 279 in the latter. 
This may imply that in 3C 279, the inverse- Compton 
(IC) 7-ray emission is generated by a slightly higher en- 
ergy population of electrons than the synchrotron opti- 
cal emission. In PKS 2155-304, the smallest timescales 
over which the optical-IR flux changed by a factor of 2 
is ~30 days, again suggesting slower variability in these 
wave bands similar to the steeper power spectral slope in 
optical-IR than that in the 7-rays as found in §3. This 
is discussed in mo r e det ails in §6. 

iTavecchio et al.l (|2010l ) showed that during a very high 
state in 2009 December, the 7-ray flux from 3C 454.3 
changed by a factor of ^5 at 6-12 hr timescale. Our 
observations of this object during that interval, as shown 
in Figure [5j had a sampling rate of only 1 — 2 times 
per night. We have obtained observations of 3C 454.3 
and PKS 1510-089 to search for intra-night variability 
during other nights. We show these light curves in Figure 
O They show significant variability over the few hour 
timescales but the amplitude of variability (~5-10%) is 
not at the same level as the hour-scale 7-ray variability 
of 3C 454.3 and PKS 1510-089 shown in ITavecchio et all 
IpOlTi) . Theintr a-night optical-IR variability shown here 
is not simultaneous with the very high state of 7-ray 
emission when 7-ray light curves with comparable time 
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Table 3 

Variability Timescales Present in the Multi-Frequency Light Curves of the Six Blazars. For the sources AO 
0235+164 and PKS 1510-089, and the 7-ray band, all timescales are upper limits. 



T,., 



(days) 



Ti, 



(days) 



Object 


B 


V 


R 


J 


K 


7-ray 


B 


V 


R 


J 


K 


7-ray 


AO 0235+164 


<1.9 


<1.9 


<1.9 


<1.9 




<1.0 


<2.0 


<2.0 


<2.0 


<2.0 




<1.0 


3C 273 












<1.0 












<1.0 


3C 279 


4.1 


4.1 


4.2 


4.2 


4.1 


<1.0 


4.9 


4.9 


4.9 


4.9 


4.9 


<1.0 


PKS 1510-089 


<2.0 


<2.0 


<2.0 


<2.0 


<2.0 


<1.0 


<1.9 


<1.9 


<1.9 


<1.9 


<1.9 


<1.0 


PKS 2155-304 


32.0 


36.0 


32.0 


33.0 


30.0 


<1.0 


23.0 


22.0 


25.0 


23.0 


25.0 


<1.0 


3C 454.3 


3.9 


3.9 


2.9 


2.9 


1.9 


<1.0 


10.0 


9.1 


12.0 


9.1 


10.0 


<1.0 



CO 2e-11 
CD 

1e-10 



5150 5155 5160 5165 5170 5175 5180 

MJD-50000 

Figure 5. Variation of 3C 454.3 at B, R, and K band during the 
large outburst in 2009 December. It is clear that the variability 
amplitude is larger at the longer wavelengths. This may be due to 
the "dilution" of the variability due to the presence of a bluer emis- 
sion from an accretion disk-like component which stays constant 
during this interval. 

resolution can be extracted with enough signal to noise 
ratio. Stronger constraints on the relation between 7- 
ray/optical-IR emission could be drawn from such data 
if they were exactly simultaneous. 

5. FLARE ANALYSIS 

Blazar light curves at all wave bands have been inter- 
preted as a superposition of outbursts caused by events 
in the jet or the accretion disk/cor ona region in addition 
to a steady baseline flux (e.g.. IValtaoia "eTaLl H99l 
i Chatteriee et alj 120081: Uorstad et alj 120101: lAbdo et al l 
2010b). Therefore, comparison of the properties of 
individual contemporaneous flares at the 7-ray and 
optical-IR wavebands is a potential dia gnostic of their 
origin . To investigate this, we follow Valtaoia et al. 
(1999) to decompose the light curves into individual 
(sometimes overlapping) flares, each with exponential 
rise and decay of the form: 

f(t) = f + f max exp[(t - t )/T r ], for t < t , and (1) 

= fit) + fma X CX P [-(t - t )/T d ] for t > t 

In this equation, f is the background level of flux that 



stays constant over the corresponding interval, f max is 
the amplitude of the flare, to is the epoch of the peak, 
and T r and Td are the rise and decay time scales, respec- 
tively. These are the five free parameters for each flare, 
/o was constrained to be less than or equal to the lowest 
value of the flux during the respective interval. There 
was no prior restriction on the possible values of the 
other four free parameters. Exponential rise and decay 
is a widely used and successful model for observationally 
well- sampled flux outbursts in blazars, at all w avelengths 
(e.g.. IValtaoia et aljfl999l: lAbdo et al.ll2010b| ). We have 
tried other models such as linear rise and decay and a 
Gaussian profile but the fit is better with the exponen- 
tial profile described in the paper. 

Similar to §2 and §3, we select the longest segments of 
R-band and 7-ray light curves of the six blazars without 
seasonal gaps (~3-4 months) to increase the accuracy of 
the flare decomposition analysis. Before the decomposi- 
tion, we smooth the light curve using a Gaussian func- 
tion with a 5-day FWHM smoothing time. Our goal is to 
compare the properties of the longer-term (longer than 5 
days) flares present in the optical light curves. We have 
shown that the PSDs of the optical variability correspond 
to red noise, i.e., there is higher amplitude variability on 
longer than on shorter timescales. Because of this, we 
analyze the more powerful longer timescale flares than 
the relatively weak flares on small time scales. The na- 
ture of the following results does not change if we use a 
3-day or a 7-day FWHM for the smoothing function. 

We proceed by first fitting the highest peak in the 
smoothed light curve to an exponential rise and decay, 
and then subtracting the flare thus fit from the light 
curve. We do the same to the "reduced" light curve, 
i.e., we fit the next highest peak. This reduces confusion 
created by a flare already rising before the decay of the 
previous flare is complete. We fit the entire light curve 
in this manner with a number of individual (sometimes 
overlapping) flares, leaving a residual flux much lower 
than the original flux at all epochs. In each case, we use 
the minimum number of flares required to adequately 
model the light curves such that using more flares does 
not change the residual flux by more than 10%. Fig- 
ures [7] and [8] compare the smoothed optical and 7-ray 
light curves with the summed flux (sum of contributions 
from all the model flares at all epochs). The 7-ray light 
curve of the blazar PKS 2155-304 was not modeled in 
the above analysis. That is because the Fermi 7-ray light 
curve of this blazar is dominated by sporadic short-term 
flares with frequent non-detection at intermediate times 
which is not suitable for fitting with the function that 
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Figure 6. The right panels show the very short-term (few hours) variability of the blazars PKS 1510-089 (J-band) and 3C 454.3 (V-band). 
The corresponding panels on the left show the longer-term variability of the same object at the same wave band. The dashed vertical lines 
in the left panels denote the intervals shown in the right panel. There is significant variability at timescales of a few hours in both objects 
but the amplitude is ~ 5 — 10%, much smaller than the 7-ray variability at these timescales. 



Table 4 

The best-fit parameters of the major individual flares in R-band. 



Table 5 

The best-fit parameters of the major individual flares in 7-rays. 



Flare Parameters 



Flare Parameters 



Object 


Flare Index 


f a 


to b 


T r c 


T d c 


e 1 


Object 


Flare Index 


f a 

<-max 


to b 


T r c 


T d c 




AO 0235+164 


opl 


4.5 


4734 


12.5 


9.0 


-0.2 


AO 0235+164 


gl 


11.0 


4767 


16.5 


20.0 


0.1 




op2 


3.9 


4753 


4.0 


20.0 


0.7 




g2 


11.0 


4722 


20.0 


20.0 


0.0 


3C 273 


opl 


28.0 


5214 


40.0 


40.0 


0.0 


3C 273 


gl 


49.0 


5094 


9.0 


17.5 


0.3 




op2 


27.6 


5270 


11.0 


8.5 


-0.1 




g2 


25.5 


5066 


12.5 


5.0 


-0.4 


3C 279 


opl 


2.1 


5040 


10.0 


17.0 


0.3 




g3 


24.5 


5144 


6.0 


6.0 


0.0 




op2 


1.1 


4949 


39.5 


31.5 


-0.1 


3C 279 


gl 


16.0 


4801 


20.0 


20.0 


0.0 


PKS 1510-089 


opl 


7.8 


4962 


7.0 


6.5 


0.0 




g2 


15.8 


5047 


20.0 


20.0 


0.0 




op2 


4.8 


4934 


16.5 


7.5 


-0.4 




g3 


14.8 


4884 


20.0 


19.5 


0.0 


PKS 2155-304 


opl 


36.0 


4966 


30.0 


30.0 


0.0 




gl 


11.0 


5023 


15.0 


4.0 


-0.6 




op2 


22.0 


5056 


17.0 


30.0 


0.3 




go 


9.7 


5078 


7.0 


20.0 


0.5 


3C 454.3 


opl 


10.2 


5170 


13.5 


19.5 


0.2 


PKS 1510-089 


gl 


35.0 


4916 


19.0 


12.5 


-0.2 




op2 


7.3 


5067 


11.0 


26.5 


0.4 




g2 


34.0 


4947 


3.5 


20.0 


0.7 
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4845 
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we use. Among all the flares used to adequately fit the 
light curves, we show the best-fit parameters of the ma- 
jor flares in Tables @] and [3] We define the major flares as 



a Flare amplitude in units of 10 7 ph cm 2 s 
b Date of peak in units of MJD-50000. 
c Rise and decay timescales in units of days, 
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Figure 7. Decomposition of continuous optical light curves of six blazars into individual flares. Red points denote the smoothed R-band 
light curves of AO 0235+164, 3C 273, 3C 279, PKS 1510-089, PKS 2155-304, and 3C 454.3. Blue solid curves correspond to summed flux 
after modeling the light curve as a superposition of several individual flares and a constant background level, magenta dotted lines show 
the residual fluxes in each case, and gray dashed lines denote the underlying individual flares. 



the ones whose amplitudes are at least at the 50% level of 
the highest peak of that light curve. We consider the the 
major flares only because in the case of the smaller flares 
the rise and decay timescales are not well-constrained. 
We define a skewness parameter (£) as the following: 



Td — T r 
T d + T r 



(2) 



£ is for exactly symmetric flares. If the decay is slower 
than the rise then £ is positive and vice versa. We show 
the distribution of the skewness parameter for all major 



optical and 7-ray flares in Figure|Hl This indicates that in 
both wave bands most of the major flares are symmetric 
(|£| < 0.3). The rise and decay times were two separate 
free parameters with no constraints. Hence, the best- 
fit models were free to contain widely asymmetric values 
of the rise and decay timescales. We showed that the 
majority of the large and well-defined flares were roughly 
symmetric. This result is not due to any constrain in our 
modeling process but comes out naturally from the data. 

The variability in the emission from blazar jets is 
primarily caused by the interaction of a shock wave 
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Figure 8. Decomposition of continuous 7-ray light curves of six blazars into individual flares. Red points denote the smoothed 0.1-300 
GeV light curves of AO 0235+164, 3C 273, 3C 279, PKS 1510-089, and 3C 454.3. Blue solid curves correspond to summed flux after 
modeling the light curve as a superposition of several individual flares and a constant background level, magenta dotted lines show the 
residual fluxes in each case, and gray dashed lines denote the underlying individual flares. Similar analysis for PKS 2155-304 was not 
carried out since in that source the variability is dominated by sporadic short-term flares with frequent non-detection at intermediate times 
which is not suitable for fitting with the function that we use. 

nite duration of the event (At) will be smoothed out 
/e.g.. iChiaberee fc Ghiselhnj[l99l iKataoka et all llM 
iSokolov et all 12004 : iChen et al.l 120111 ). If the radiative 
cooling time is longer than At, then we would expect to 
see a significantly longer decay time than rise time. This 
is because the radiative cooling timescale of electrons 
which are generating the optical-IR emission through 
synchrotron radiation and 7-rays through IC processes 
is longer than time needed to inject particles (usually 
assumed to be "instantaneous" ) responsible for the flare 
into the region where they would be able to radiate sig- 



with the jet plasma (e.g., iBlandford fe Konigl 119791 : 
IMarscher fc GearlfTflS^ over a finite size or a finite time 
interval. The former can occur when a disturbance in 
the jet-flow or a blob of denser plasma, observed as mov- 
ing knots in radio-interferometric observations, passes 
through a standing shock present in the jet. The lat- 
ter can occur due to the evolution of an internal shock. 
Due to the existence of a finite size and the resultant 
finite duration of the above events, the observations are 
affected by light-travel time. Any process faster than 
the light-travel time through the said size or the fi- 
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Figure 9. Distribution of the skewness parameter (£) for the ma- 
jor optical and 7-ray flares. This shows that most of the large flares 
are symmetric in nature (|£| < 0.3) at both wave bands. 

nificantly. On the other hand, we may expect to see 
symmetric flares, i.e., the rise and decay times are com- 
parable, if the radiative cooling time of the emitting par- 
ticles is much smaller than At. For example, in the case 
of a disturbance in the jet moving through a standing 
shock, the rise time corresponds to the time it takes for 
the entire disturbance to enter the standing shock re- 
gion while the number of energized electrons increases 
resulting in an increase in flux. Similarly, the decay time 
corresponds to the time taken by the disturbance to leave 
the standing shock region entirely while the number of 
energized electrons steadily decreases since the radiative 
cooling times are much smaller than the rise and decay 
times. In this case, the rise and decay times will be simi- 
lar resulting in a symmetric flare. The symmetric nature 
of the observed flares in this work indicates that the rise 
and decay timescales are dominated by crossing time of 
radiation or a disturbance through the emission region, 
similar to the latter case described above. 



6. 



DIFFERENCE IN THE TIME VARIABILITY 
PROPERTIES OF FSRQS AND BL LACS 



The main difference between the time variability prop- 
erties of these two classes of objects that has been re- 
vealed in this work is that the R-band variation is signif- 
icantly smoother than that in the 7-ray energies in the 
only BL Lac object in our sample, namely, PKS 2155- 
304 while in the FSRQs the variations in these two wave 
bands are similar. This may be caused by one or both of 
the following reasons: 

(1) The 7-ray emission in blazars is generated by the 
IC processes. Hence, the variability is a combination 
of the variation of the emitting electrons and that of 
the seed photons which are being inverse- Compton scat- 
tered to the 7-ray energies. In the case of the HBLs 
such as PKS 2155-304, GeV 7-ray emission is dominated 
by the SSC process while in the FSRQs it is dominated 
by the EC contribution. In the former, the seed pho- 



tons which are produced by the synchrotron emission in 
the jet are varying faster than those in the latter which 
are generated in the disk, broad line region or the torus. 
Hence, in the case of PKS 2155-304, combined variability 
of the emitting electrons and the seed photons which is 
reflected in the GeV band is faster than that in the R- 
band which represents the variability of the synchrotron 
emission only. This is not seen in the FSRQs where lit- 
tle or none of the GeV variability is due to the variation 
in the slowly varying seed photons. Hence the GeV and 
R-band variability is very similar because both are con- 
tributed by the variation in the emit ting electrons only . 
(2) From the SED of PKS 2155-304 (jAbdo et alJl2lTi0dY 
it can be seen that the R-band emission is ~1.5 orders 
of magnitude lower than the synchrotron peak while the 
7-ray peak is at 4 GeV, which is within the wave band 
we are considering in this paper, 0.1-300 GeV. The pho- 
ton index of this object from the Fermi 2-yr catalog 
(jAbdo et al.H20Tl is 1.84. This means in this blazar, the 
0.1-300 GeV emission has significant contribution from 
electrons which are at a higher energy than those produc- 
ing the synchrotron radiation unlike the FSRQs where 
the GeV and R-band emission is produced by electrons 
of similar energies. This may partially contribute to the 
faster variability in the GeV band than that in the R- 
band in PKS 2155-304, as observed. 

7. DETAILED COMPARISON OF THE 7-RAY AND 

OPTICAL LIGHT CURVES DURING A PROMINENT 
MULTI-FREQUENCY OUTBURST OF THE BLAZAR 3C 
454.3 

The flat spectrum radio quasar 3C 454.3 went through 
a very large multi-wavelength outburst during 2009 De- 
cember. There were several smaller but significant flares 
at 7-ray energies as well as optical-IR frequencies during 
the four months prior to the large outburst. To investi- 
gate the location of these outbursts at 7-ray energies as 
well as to compare their properties with corresponding 
flares at optical-IR frequencies, if present, we decompose 
the light curves of 3C 454.3 from Fermi and SMARTS 
into individual, sometime overlapping flares in the same 
manner as in §5. 

To ensure that even small-timescale flares are resolved, 
wc smooth the 7-ray and optical-IR light curves of the 
blazar 3C 454.3 during 2009 August-December with a 
Gaussian smoothing function of FWHM 2 days. Fig- 
ure [10] compares the smoothed light curves with the 
summed flux (sum of contributions from all the model 
flares at all epochs). We identify 6 7-ray/optical flare 
pairs in which the flux at both wavebands peaks at the 
same time within ±3 days. For each of the flares, we 
determine the time of the peak, width (defined as the 
mean of the rise and decay timescales), and area under 
the curve from the best-fit model. We calculate the area 
under the curve for each flare by integrating the photon 
flux over 0.1-300 GeV u sing average spect ral index from 
the Fermi 2-yr catalog (|Abdo et al.l 1201 If) and then in- 
tegrating over the duration of the flare. This represents 
the total energy output of the outburst. Table [6] lists the 
parameters of each flare pair, along with the ratio, p, of 
7-ray to optical energy output. 

The presence of a corresponding R-band flare for each 
of the significant 7-ray flare implies that the emission 
at both wave bands is generated by the same electrons 
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Figure 10. Red open circles show smoothed 7-ray (0.1-300 GeV) 
and optical (R-band) light curves of the blazar 3C 454.3 during the 
second half of 2009. Blue solid curves correspond to summed flux 
after modeling the light curve as a superposition of six individual 
flares and a constant background level while magenta dotted lines 
denote the residual flux after subtracting the model flares, and the 
gray dashed lines indicate the underlying individual flares. Flares 
denoted by the same number in both panels peak within ±3 days 
of each other and are probably physically related ("flare-pairs"). 
These flare pairs are listed in Table [6] 

and h ence provides strong support to the leptonic models 
(e.g.. iBonning et al.ll2009f) . However, with some tuning, 
some hadronic models can also reproduce this behavior, 
so our results do not definitively rule out this alterna- 
tive. The total energy output of the 7-ray flares are 
larger than that in the optical by 1-2 orders of magni- 
tude on average. This is con sistent with the SED of this 
source (|Bonning et a.1.1 [2012T ) which shows that it emits 
much greater energy at 7-ray-frequencies than that at 
the optical wave band. Table |H] also shows that the en- 
ergy output at the same energy band as well as the total 
energy output of the sum of a 7-ray flare and the corre- 
sponding optical flare varies significantly from one event 
to the other. Figure [Til shows the 7-ray flux vs. optical 
flux plot for each of the identified flare pairs. The lines 
denote three forms of the numerical relation between the 
fluxes in those two bands, e.g., F 7 ~ F* p , where x=l, 
1.5, and 2. 

Exact location of the 7-ray production region in blazars 
is not well-determined. Optical-UV photons coming from 
the broad line region (BLR) may be inverse- Compton- 
scattered by the electrons in the jet, in which case the 
bulk o f the observed lumino sity appears in the 7-ray band 
(e.g- IGhisellini et al.|[2010h . But previous time variabil- 
ity studies from the EGRET era have indicated that the 
7-rays are produced downstream of the very long baseline 



interfero metry (VLBI) core (which lies >1 pc from the 
SMB H) (porstad et all 120011: iLahteenmaki fc Valtaoial 
2003). 7-7 absorption of very high energy 7-ray pho- 
tons by the radiation field of the BLR may present 
another problem for a model of very high energy 7- 
ray emission inside the BLR of lumino us quasars (e.g., 
iDonea fc Protherod 12001 iReimerl l2007f ). On the other 
hand, the external photon density at >1 pc from the 
base of the jet is theoretically not large enough to gen- 
erate the amount of 7-rays observed from many blaz ars 
(e.g., iGhisellini fc Tavecchiol[2009t ISikora et al.fl2009l ). 

Recently. lAgudo et al.1 f|201 lh have shown that bright 
7-ray flares in the jet of the blazar OJ 287 occur >14pc 
from the central engine. They find this by analyzing a 
combination of time-dependent multi- waveband flux and 
linear polarization observations, and sub-milliarcsecond- 
scale polarimetric VLBI images at 7 mm. Similar lo- 
cation (>few pc from the black hole) of the production 
of 7-r ay e mission has been infer red by IMarschc r et al.l 
( 2008) and lMarscher et al.l ()2010D as well for the blazars 
BL Lac and PKS 1510-089, respectively. They argue that 
large multi-waveband outbursts are triggered by the in- 
teraction of moving plasma blobs with a standing shock 
present in the jet seen as the "core" in the pc-scale jet in 
Very Long Baseline Array (VLBA) images. In all these 
cases, the large multi-frequency flares were preceded by 
smaller flares caused by the movement of those plasma 
blobs along the jet axis. The blobs follow a helical path 
at the acceleration-collimation zone of the jet due to the 
nature of the magnetic field present there. That causes 
the movement of the blobs to align with our line of sight 
and hence their emission to be Doppler boosted during 
some intervals. 

To understand the location and mechanism of the 
7-ray and optical emission in this blazar we discuss the 
dependence of the observed synchrotron (F sync h), syn- 
chrotron self-Compton (Fssc)j and external-Compton 
(Fec) fl ux on three relevant parameters, namely, total 
number of emitting electrons (N c ), magnetic field (B) 
and Doppler factor (<5). We choose these parameters 
because these are independent of each other in the 
following functions: 

F synch ~ N B 1+ «°5 3+Q ° (3) 
F EC ~ N c <5 4+2 ^U; xt (4) 
F S sc ~ N c <5 3 +^U^ ynch , (5) 

where qq and a g are the spectral indices of the syn- 
chrotron and IC emission respectively, Ug Xt is the ex- 
ternal seed photon field, and Ug ynch is the same due to 
synchrotron emission in the jet itself. The observed syn- 
chrotron emission from movi ng plasma in a relativisti c jet 
is amplified by a factor S 3+a (|Urrv fc Padovanilll995f ) . In 
the case of EC emission the amplification factor is S 4+2a 
due to the additional Lorentz tra nsformation b etween the 
seed photon and jet rest frame (|Dermerl ll995). There is 
no additional factor of S in Fssc in Equation 5 because 
the synchrotron seed photons are also from the jet rest 
frame where the emitting particles reside. We note that 
i) N e is the number of electrons contributing to the ob- 
served flux and is related to the normalization (No) of the 
electron spectrum. The actual dependence of flux is on 
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Table 6 

Time, Total Energy Output (Area), and Temporal Width of 7-ray/Optical 
Flare Pairs Identified in the Light Curves of the Blazar 3C 454.3 During 
2009 August-December. 



Flair Pair 7-ray optical p d 



ID 


Date a 


Area b 


Width c 


Date a 


Area b 


Width c 




1 


5071 


180 


6.5 


5068 


5.6 


10.5 


32 


2 


5090 


356 


7.5 


5089 


2.1 


7.0 


170 


3 


5112 


164 


11.0 


5109 


4.7 


12.8 


35 


1 


5127 


380 


11.5 


5125 


1.0 


7.0 


380 


5 


5148 


13 


1.8 


5147 


0.2 


1.5 


65 


6 


5169 


1084 


10.5 


5168 


6.5 


18.3 


167 



a Date of peak in units of MJD-50000. 
b Units: 10 -5 erg cm -2 . 

c Mean of rise and decay timescale in units of days. 

d Ratio of 7-ray to optical energy output integrated over flare. 



No. If the minimum energy of the electron distribution 
(7min) remains unchanged N e is proportional to No, and 
ii) Above mentioned fluxes are at fixed observed energy 
bands as shown in Figure [TT] There are other relevant 
quantities such as the size of the emitting region but we 
assume that it remains constant in order to determine the 
effect of the above parameters on the observed flux. A 
detailed numerical calculation is required to accurately 
model the effect of various quantities on the observed 
fluxes but such a calculation is b e yond the scop e of the 
paper (see iTramacere et al.1 120091 : iDermerl 119951 . for de- 
tails). Our goal here is to compare the observed data to 
approximate theoretical scenarios. 
From Equation 3 and 5, we get: 

F ssc ~ N lB 1+ao 5 3+a z (6) 

Therefore, i) If the variation is due to a change in B, 
Fec will not vary while F syncn ~B 1+a ° , i.e., the 7-ray 
and the optical variation will not be correlated. On the 
other hand, Fssc^Fsynch h) If it is due to a variation in 
N e , F E c~F S ynch and F S sc~Fsynch and in ) If tne variation 
is due to a change in S, F EC ~ F^+ 2 h Qs)/(3+Qo) and F SSC ~ 

Fsynch 9 ^ 3+Q °^ During the interval of these flares, i.e., 
MJD 55060 to 55160, OIR spectral index from SMARTS 
light curves is 1.55±0.05. The 7-ray spectral index dur - 
ing the same interval is 1.5±0.1 (|Ackermann et ai] [2010). 
For a = 1-55 and a g = 1.5, F E c ~ F^ ch and F S sc ~ 
F^ynch- Therefore, from Figure [Til we can confirm that 
if the 7-rays are EC in nature, the variation is not due 
to a change in B. No other possibilities can be ruled out. 
Other relations between the 7-ray and optical flux vari- 
ation is possible if the variation is due to a combination 
of two or more of the above scenarios. 

It is possible that the large multi-frequency flare of 
3C 454.3 during 2009 December is another case where 
a large outburst is generated due to the interaction of a 
moving plasma knot with a standing shock present in the 
jet. This can be tested by following the variation in the 
structure of the pc-scale jet with VLBA imaging for a suf- 
ficient duration. Such a studj0 indeed shows that a new 
knot (moving plasma blob) was coincident with the core 

6 http:/ /www. bu.edu/blazars/VLBAproject. html 
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Figure 11. The data points show the 7-ray vs. optical fluxes for 
each of the flares identified in Tab le [6 ] with each symbol denoting 
a different flare number in Figure [T0l magenta asterisks-1, green 
filled circles-2, cyan triangles-3, blue filled squares-4, yellow filled 
circles-5, and red open circles-6. The lines denote three forms of 
the numerical relation between the fluxes in those two bands, e.g., 
F 7 ~ FJ p , where x=l (magenta dotted), 1.5 (blue dashed), and 2 
(red solid). 

(standing shock) on 2009 November 27 (±15 Days) (Alan 
Marscher, private communication). This implies that the 
large 7-ray and optical outbursts might take place in the 
jet near the VLBA core located at ~18 pc from the cen- 
tral engine ([Jorstad et a,l.M2010h . 

In that case, the smaller flares preceding the large out- 
burst may be due to Doppler boosting of the emission 
from the same knot when it is propagating through the 
acccleration-collimation zone of the jet and getting in and 
out of our line of sight. This is consistent with the fol- 
lowing properties of these flares: i) We can see in Figure 
llll that the 7-ray vs. optical flux data points for four out 
of five smaller flares (yellow, green, cyan, and magenta) 
are roughly consistent with F 7 ~ Fj p 5 or F 7 ~ F^-p, 
as expected if the variation is due to changes in the 
Doppler factor, ii) The range of F 7 values on Figure [TT] 
corresponding to the five flares preceding the large flare 
around MJD 55169 denoted by flare number 6 in Fig- 
ure [10] correspond to a change by a factor of ~3.7 while 
that for F op is ~2.8. The fractional change required in 5 
to accoun t for these changes i s 20-25% which is definitely 
plausible ([Jorstad e t al. 2005) and finally iii) These flares 
were near-simultaneous at all wave bands and the recur- 
rence time of the peaks were similar. The first peak is 
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on MJD 55071 and after that each flare is ~20 days after 
the previous flare within ±3 days. This can be explained 
by the movement of the plasma blob through a helical 
path while its emiss ion is affected by Doppler boosting 
at regular intervals (iCamenzind fc Krockenbergerlll992l : 
iMarscher etl!ll2lM 120101 : lAgudo et alJl2Qlll ). 

The scenario described above is not the only possible 
process to explain the above observations. As described 
above, the variation could be due to changes in N e , S 
or a combination of one or both of them with changes 
in B. The physical reason for a change in N e or B is 
not clear from the data. For example, the acceleration 
of electrons due to the presence of instability in the jet 
could increase N e . Alternatively, stochastic injection of 
energy in the jet plasma from the central engine and its 
dissipation at differe nt distances from the base of the jet 
(Bonn oli et al.l 120111) could explain the observations as 
well. 

If the Doppler beaming explanation is true, this adds 
to the increasing evidence that at least in some blazars, 
copious 7-ray emission is produced farther down the 
jet and in those cases the BLR cannot be a significant 
source of seed photons which are inverse- Compton scat- 
tere d to 7-ray energ i es by the energetic electrons in the 
jet. iTavecchio et al l ((20100 showed that during this high 
state in 2009 December, the 7-ray flux from 3C 454.3 
changed by a factor of ^5 on a 6-12 hr timescale. This 
implies that the size-scale of emission region, obtained 
from the variability timescale (~0.01 pc) is smaller than 
the approximate cross-section of the jet o f 3C 454.3 at 
18 pc (~0.05 pc, assuming 0j Ct ^O.2 o from Uorstad et afl 
2005). This may imply that the size-scale of the 7-ray 
emission regions is smaller than the diameter of the lo- 
cal cross-section of the jet. Both of these provide strong 
constraints to the theoretical models of 7-ray production 
in the jets of blazars. 

Stronger constraints on the location and mechanism of 
emission will come from a detailed SED modeling of the 
SMARTS, Fermi and other multi-wave band data of the 
above individual successive flares, to be addressed in a 
future paper. 

8. SUMMARY AND CONCLUSIONS 

This paper presents the time variability properties of a 
sample of six 7-ray-bright blazars at optical-IR frequen- 
cies as well as 7-ray energies: AO 0235+164, 3C 273, 3C 
279, PKS 1510-089, PKS 2155-304, and 3C 454.3. The 
light curves were obtained as part of the Yale/SMARTS 
program in 2008-2010 to monitor all bright southern 
Fermi-LAT-monitored blazars on a regular cadence, at 
optical and near-infrared (BVRJK) wave bands. Our 
main conclusions are as follows: 

(1) We find the optical-IR variability properties to be 
remarkably similar to those at 7-ray energies. This is 
consistent with the general picture of the leptonic model 
where the lower (optical-IR) and higher (7-ray) energy 
emission is generated by the same population of electrons 
through synchrotron and inverse- Compton processes, re- 
spectively. However, more rapid variations indicate that 
the electrons producing 7-ray emission in the blazar 3C 
279 may have slightly larger energy than those gener- 
ating optical-IR radiation in this source. We note that 
some hadronic models are also able to reproduce the sim- 
ilarity of the optical-IR variability properties of blazars 



with those at the 7-ray bands and hence we cannot rule 
out those models on this basis. 

(2) The discrete auto-correlation function (ACF) of the 
variability of these six blazars at optical-IR and 7-ray 
wave bands do not show any periodicity or characteristic 
timescale. The shape and width of the ACFs are very 
similar in all bands in the 3 sources where the temporal 
sampling at both energies are identical. This indicates 
that the emission regions are of similar sizes with light- 
crossing time being the dominant timescale. 

(3) The power spectral density functions of the R-band 
light curves of all six blazars are fit well by simple power- 
law functions with negative slopes, implying there is 
higher amplitude variability on longer than on shorter 
timescales. No evidence of a characteristic timescale, 
including a break, is identified in the PSD of any of 
the sources. The slope of the aver age PSD of the sam - 
ple of six blazars is similar to what lAbdo et al.l ()2010bl ) 
found for the 7-ray variability of a larger sample of bright 
blazars. This is consistent with the conclusion in item 
(1) that the lower (optical-IR) and higher (7-ray) energy 
emission is generated by the same population of elec- 
trons. 

(4) For all six sources, the shortest timescales over which 
the 7-ray flux changes by a factor of 2 is <1 day. For the 
blazars AO 0235+164 and PKS 1510-089, this timescale 
in optical-IR energies is 1 — 2 days, consistent with those 
in the 7-ray band. However, in 3C 454.3, the dou- 
bling timescale is longer at shorter wavelength optical- 
IR bands. This may be due to the disk contribution 
in the bluer bands which dilutes the jet variability and 
hence increases the doubling timescale. In 3C 273 the 
flux does not change by a factor of 2 during this observ- 
ing program which may also be a result of significant 
contribution from the disk in the optical-near IR emis- 
sion of this object. 

(5) The variability properties, including the ACF peak 
value, the relative nature of the short and long-term vari- 
ability at 7-ray and optical wave bands, and the doubling 
timescales are different in the BL Lac 2155-304 from the 
FSRQs in our sample. This may be caused by: i) the 
7-rays in PKS 2155-304 are dominated by the SSC pro- 
cess unlike the FSRQs in which it is dominated by the 
EC contribution and ii) the electrons generating the 7- 
ray emission in this blazar are at a higher energy than 
those producing the optical-IR emission unlike the FS- 
RQs where the GeV and R-band emission is produced by 
electrons of similar energies. 

(6) We decompose the optical-IR and 7-ray light curves 
to investigate the properties of the individual flares 
present in those light curves. The prominent flares 
present in the optical-IR as well as the 7-ray light curves 
of these blazars are predominantly symmetric, i.e., have 
similar rise and decay timescales. This indicates that the 
the long-term variability is dominated by the crossing 
time of radiation or a disturbance through the emission 
region and not by the energy-loss time of the emitting 
electrons due to radiation losses. 

(7) Finally, we perform a detailed decomposition of the 
SMARTS and Fermi light curves of 3C 454.3 during and 
before its large outburst in 2009 December to understand 
how and where the emission in the prominent outbursts 
of blazars is generated. The total energy output, the ra- 
tio of 7-ray to optical energy output, and the 7-ray vs. 



15 



optical flux relation of six individual flares of the blazar 
3C 454.3 during 2009 August to December vary signifi- 
cantly from one event to the other. 

(8) In 3C 454.3, the smaller flares at both optical-IR and 
7-ray wave bands during 2009 August-November, and 
the subsequent large multi-frequency outburst in 2009 
December, might be triggered by the movement of a 
plasma blob along a helical path in the jet and its sub- 
sequent interaction with a standing shock present in the 
jet, respectively. If this scenario is right, the location 
of the large 7-ray outburst in 2009 December is iden- 
tified at ~18 pc down the jet from the central engine, 
where a strong radio-emitting component is observed. 
This would pose strong constraints on the models of high 
energy emission in the jets of blazars. 
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